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Abstract 

In this work, the edge physics of an Aharonov-Bohm interferometer (ABI) defined on a two dimensional electron gas, subject 
to strong perpendicular magnetic field B, is investigated. We solve the three dimensional Poisson equation using numerical 
techniques starting from the crystal growth parameters and surface image of the sample. The potential profiles of etched and 
gate defined geometries are compared and it is found that the etching yields a steeper landscape. The spatial distribution 
of the incompressible strips is investigated as a function of the gate voltage and applied magnetic field, where the imposed 
current is confined to. AB interference is investigated due to scattering processes between two incompressible "edge-states". 
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The recent increasing interest towards the quantum 
Hall based interferometers relies on the popularity of 
the quantum information processing. In particular, the 
realization of electron and quasi-particle interference 
experiments became a paradigm [1,2]. A possible ap- 
plication is proposed to use the non-Abelian 5/2 state 
for topological quantum computation [3], which es- 
sentially has a very similar structure of AB interfer- 
ometers. The well established experimental wisdom is 
that, to realize "clean" measurements are extremely 
difficult, which strongly depends on the sample geom- 
etry, crystal growth etc. In typical AB interference ex- 
periments two propagating states are brought to close 
vicinity, by the help of gates [1] . The edge states form 
a closed (or almost closed) path, which enclosures cer- 
tain amount of magnetic flux. By changing the mag- 
netic field or the area of the closed path, one infers 



the phase of the particles. The conventional edge pic- 
ture is used to explain the observed AB oscillations 
[4], however, the actual distribution of the edge-states 
is still not known for realistic samples, although, sev- 
eral powerful techniques are used [5]. At the recent 
experiments of Camino et. al [2], show that the con- 
ventional AB theories, which neglect electron-electron 
interactions, are unable to explain the periodicity of 
the oscillations. The only very recent model of Igor 
Zozoulenko [5] could provide a reasonable explanation 
for the unexpected behavior, which states that "stan- 
dard transport models" fail to understand the underly- 
ing physics. However, the model geometry considered 
in their work is quite different from the actual experi- 
mental setup. Here, we provide an explicit calculation 
scheme to obtain the density and potential profiles of 
an AB interferometer in the absence of an external 
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Fig. 1. (a) and (c) considers only the gated sample, whereas 
b) and d) is for both etched and trench gated structure, 
a) and b) simulate the potential profiles of the sample and 
the corresponding electron densities are shown in (c) and 
(d). We ste the same gate voltage (-1.8 V, note that surface 
potential is -0.7 V due to mid-gap pinning) for both situ- 
ations. Potential distribution of the gated sample presents 
smooth characteristic, however, the potential profile of the 
trench gated structure is very sharp. Hence the electron 
density at the island center is higher in (d) compared to 
(c), in the island center. All the figures are obtained at zero 
temperature and magnetic field. 

magnetic field and also under quantized Hall condi- 
tions. Our calculation is based on solving the Poisson 
equation in 3D, starting from the crystal growth pa- 
rameters and the lithographically defined surface pat- 
terns. A fourth order nearest neighbor approximation 
is used on a square grid with open boundary conditions 
and 3D fast Fourier transformation method is used to 
obtain the solution iteratively [6,7]. The outcome, i.e. 
the potential and electron density distributions, of this 
calculation is used as an initial condition for the mag- 
netic field dependent calculations. The distinguishing 
part of our calculation is that we do not have to as- 
sume only gate defined structures but we can also han- 
dle etching defined geometries, which essentially is the 
case for the experiments. We show that, the etching 
defined samples present a sharper potential profile [7] , 
therefore, the formation of the edge-states are strongly 
influenced. This, obviously, effects the edge physics in 
determining the AB oscillations. 

The simple description of an ABI, is such: let us 
assume a closed path, a circle for simplicity, subject to a 
perpendicular magnetic field with an area of yrr^ where 
r is the radius. Since the wave function travelling from 
one side should be the same as the one the travelling 
from the other side, the phase difference of the wave 
functions can be only integer multiples of the magnetic 
flux quanta $0 = h/e encircled. In other words, for a 



given flux $ every spin resolved cyclotron orbit should 
satisfy the condition $ — BSm ~ rn^o, where m is 
the quantum number of the orbit. Therefore, an orbit 
with a radius of magnetic length I = {^/n/eB ) will 
have an area Sm = 27vml^, resulting in Sm+i — Sm ~ 
h/eB. It is common to define the occupation of these 
orbits by = n^h/eB, where n^i is the number of 
the electrons and this occupation is called the filling 
factor. This picture should also hold for non- interacting 
electrons which are confined by an external potential 
14xt(r), which essentially lifts the degeneracy. In the 
experiments considered in this work, it was shown that 
the number of electrons is fixed (nd = 1700 — 2000, 
varying with the sample size) in the "island", which 
means that it is energetically almost impossible to add 
or subtract an electron from the quantum dot region. If 
one keeps the electron density fixed and decreases the 
B field by a factor of two, the filling factor increases by 
a factor of two, which implies that the are enclosed by 
= 1 is now half. 

In this work we only consider the ABI experiments 
conducted at the integer quantize Hall regime, i.e. 
v = Uai/nB ~ k and k = 1,2,3,... [1] and the re- 
sults reported in Ref. [2] at the integer regime. The 
first work concentrates the AB oscillations observed 
only at the v — 1 plateau. In the above mentioned 
experiments, the magnetic field B and the front gate 
VpG dependency of the AB period is investigated and 
the results are discussed in terms of the Gelfand and 
Halperin (GH) [8] and Chklovskii, Shklovskii and 
Glazman (CSG) [4] edge models. A hybrid formula 
is then used to describe the actual electron density 
distribution non- self- consistently. Procedure is as fol- 
lows: The GH model describes (almost) properly the 
etched edge density profile and CSG model provides a 
Vfg dependency. Therefore, the density distribution 
without gates is taken from the GH model and its evo- 
lution depending on the front gate bias is described by 
the CSG model. It is argued that, this hybrid model is 
in agreement with the experiments with a difference of 
13% when comparing the surface area change AS^ of a 
single edge-channel at the gap as a function of Vfg- 
Such a relatively small difference, at a first glance, 
looks impressive. However, in the second report [2], 
where the results at = 2 is also shown, it was stated 
that if the radius of the outer ring remains unchanged 
the AB oscillations can be explained. 

Our investigation is based on the calculation of 



2 




"0 12345012345 
X (urn) X ((im) 



Fig. 2. Spatial distribution of the incompressible strips [(a) 
and (b)] and current distribution [(c) and (d)] are shown 
at S = 7.8 T [(a) and (c)] and 8.6 T [(b) and (d)]. Blue 
highlight the ISs (u = 2) and shaded areas indicate the 
depleted regions in (a) and (b). Blue arrows indicate the 
current distribution in (c) and (d). Fixed external current 
is driven in x-direction where the ID current density is set 
to be —1.06 X 10~* A/m. Here we only consider the trench 
gated structure (140 nm etched and VpQ = —1.8 V). 

the electron density and screened potential self con- 
sistently within the Thomas-Fermi approximation 
(TFA) [7]. We consider spinless electrons in the high 
magnetic field regime, consequently the effective 
Hamiltonian reads, 

H^^'^' = Ho + Kxt(x,y) + Vint(x,y). (1) 

Here Ho is the kinetic part, Vcxt{x,y) and Vint{x,y) 
are the external and interaction potentials, respec- 
tively. The external potential is obtained by the above 
mentioned 3D calculations considering the real exper- 
imental structure, whereas the interaction potential 
(Hartree potential) is calculated for a given density, 
boundary conditions and gate pattern by solving the 
Poisson equation 

Vint{x,y)^— / K{x,y,x',y')nci{x',y')dx'dy' (2) 

where k is the dielectric constant (= 12.4 for GaAs). 
The Kernel K{x, y, x' , y') is the solution of the Poisson 
equation preserving the periodic boundary conditions. 
Spatial distribution of the electron density is obtained 
within the TFA via. 



n^iix, y) = J dED{E)f{E + V{x, y) - ^*), (3) 

where D{E) is the relevant (collision broadened) den- 
sity of states, f{E, ^* , T) determines the particle statis- 
tics (Fermi function) and ^* is the electrochemical 
potential. The total potential V{x,y) = Vin%{x,y) -\- 
Vcyit{x,y) and eq. 3 are calculated self-consistently to 
obtain the electron and potential distributions at finite 
magnetic field and temperature. We first compare the 
results of the gate and trench gate defined samples. We 
show the potential profiles and electron distributions 
both of these samples in Fig. 1, the upper panels illus- 
trate the potential profiles and it is seen that the trench 
gate defined sample has a sharper profile. Fig Ic and 
fig. Id represent the electron densities with gray scale, 
calculated for two different definitions. These results 
are obtained at zero magnetic field and temperature 
as an initial condition for our non-zero magnetic field 
calculation. The maxima of the electron density at the 
center of the q-dot matches perfectly (< %1) with the 
experimental value when considering the trench gated 
numerical simulation. 

Next we calculate the potential and electron den- 
sity distribution at finite temperature and B field, only 
considering the trench gated sample. In Fig. 2 we show 
the spatial distribution of the ISs at two different mag- 
netic field values. Shaded areas indicate that the elec- 
trons on the 2DEG are depleted perfectly. The side- 
surface electrons which exists on the side surface due 
to etching, yield larger shaded regions than the only 
gated model. Since the potential is sharper at trench- 
gated defined samples the widths of the incompressible 
strips are narrower compared to gated samples. When 
considering the B field regime 2 < v < A there is only 
a narrow interval where two incompressible strips are 
close to each other fig. 2a, however, do NOT merge 
2b. Fig. 2c and fig. 2d presents the corresponding cal- 
culated current distribution utilizing the local Ohm's 
law [9,10,7]. We clearly observe that, the imposed fixed 
current is confined within the incompressible regions, 
where backscattering is absent. It is known that the 
current flowing from the incompressible strips is di- 
vergent free, thus it is not possible to inject current 
directly to these strips. Hence, without scattering be- 
tween two "edge channels" one would not be able to 
observe any interference pattern. The scattering mech- 
anism comes from the impurities or the electric field 
at the boundary between the ISs and compressible re- 
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gion [11], which we implicitly include to our calcula- 
tions when calculating the local conductivity tensor 
elements using the findings of self-consistent Born ap- 
proximation [10]. 

We aim to investigate a situation where, the spin- 
degenerate groun state Landau wave functions (LWs) 
overlap at a certain spatial region where scattering ma- 
trix elements become finite, hence we can observe the 
Aharonov-Bohm interference (ABI). Since the screen- 
ing is merely poor within the ISs, the total potential 
exhibits a variation where the electron density is con- 
stant [12]. Therefore an electric field develops within 
this regions if we have a slope on the potential profile, 
that can be computed from Ex = Using this 

equation and considering the calculated total potential 
profile on can obtain the spatial shift, Xe, of the center 
coordinates Xo of the LWs from 



note that the LWs remain unaffected (i.e a Gaussian 
centered at Xq) in the compressible regions, since the 
potential is (almost) constant. Now, if the ISs are large 
and sufficiently far apart, no scattering processes can 
take place hence no interference can be observed. Such 
a case is shown in fig. 3a. Whereas, if the two incom- 
pressible edge states are close enough to each other and 
Ex is large and LWs overlap, i.e. scattering takes place. 
This is shown in Fig. 3b. If we apply a high magnetic 
field ISs become broaden contrarily Ex take small val- 
ues so the shift of the center coordinate is negligible 
and as a result no overlap occurs. Since the widths of 
the incompressible strips are related to the pattern ge- 
ometry and crystal growth parameters, the observation 
of AB interference patterns are extremely fragile. For 
such experiments, one should design the sample geom- 
etry and choose the magnetic field interval keeping in 
mind that the formation and the spatial distribution 
of the incompressible strips are important. The calcu- 
lation of the actual scattering matrix elements for the 
real structures considering self-consistently calculated 
wave functions is beyond the scope of the present pa- 
per. However, we have drawn the calculation scheme. 

To summarize: We have calculated self-consistently 
the electron and potential distribution considering the 
experimental sample geometry and material proper- 
ties. It is shown from the 3D calculations that, the 
trench gated structure can be simulated in a fairly good 
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Fig. 3. The (screened) potential and the shift of the LWs 
for different magnetic fields. Slope of the potential is small 
at higher magnetic field (a) compared to low B value. 

agreement (with %1 error). In the second part of our 
work we have also considered the interference condi- 
tions depending on the magnetic field value and ob- 
served that the scattering from one edge to the other is 
only possible if the two incompressible strips are close 
enough, i.e. at the order of magnetic length. If the two 
"edge states" merge no partition takes place, hence in- 
terference pattern is smeared. 
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